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Abstract. We discuss the subject of Unextendible Product Bases with the or- 
thogonaUty condition dropped and we prove that the lowest rank non-separable 
positive-partial-transpose states, i.e. states of rank 4 in 3 X 3 systems are al- 
ways locally equivalent to a projection onto the orthogonal complement of a 
linear subspace spanned by an orthogonal Unextendible Product Basis. The 
product vectors in the kernels of the states belong to a non-zero measure sub- 
set of all general Unextendible Product Bases, nevertheless they can always 
be locally transformed to the orthogonal form. This fully confirms surprising 
numerical results recently reported by Leinaas et al. Parts of the paper rely 
heavily on the use of Bezout's Theorem from algebraic geometry. 



1. Introduction 

Unextendible Product Bases (UPBs for short) consisting of orthogonal vectors 
are well recognized by the quantum information community [IHZ] and among mathe- 
maticians [8]|^ , but similar product bases with the orthogonality condition dropped 
seem to have received very little attention [IOHl2] so far. This is partly explained by 
the fact that an orthogonal UPB can always be used to construct a PPT bound en- 
tangled state [T], while no simple prescription is known for general UPBs, cf. [TU1|13| . 
Nevertheless, the recent results of [TT], if not plain curiosity, seem to call for a more 
systematic treatment of the subject of general UPBs. It has been observed in numer- 
ical searches [TT] that non-separable positive-partial-transpose states (PPT states 
for short) of rank four in 3 x 3 systems can be locally, i.e. by SL (3, (D) (g) SL (3, (D) 
tranformations, brought to the form of an orthogonal projection onto the orthogo- 
nal complement of an orthogonal 3x3 UPB ■ The product vectors in the kernel 
of such PPT states belong to a non-zero measure subset of the set of all general 
3x3 UPBs with a minimal number of elements. A valuable discussion of the subject 
has been included in the recent works [12l[l3]. In the present paper, we provide an 
analytical proof of the main result of and explain a number of other properties 
of extreme PPT states, obtained numerically and reported in [T3]. As a byprod- 
uct, we deliver a number of rather general techniques to tackle problems related to 
product vectors in the kernel of a PPT state. Facts from algebraic geometry, and 
specifically the Bezout's Theorem are of much use here. A recognized textbook on 
the subject is [15], and we are going to use it as a main reference. 

The content of the paper is divided into two sections, each of them consisting 
of smaller pieces. The first section discusses the subject of general Unextendible 
Product Bases, without any constraints on dimensionality, neither on the number 
of subsystems. It also summarizes what is generally known for 2 x n systems 
as a consequence of the recent work [16j . and partly also by the early findings 
of [1]. The second part of the paper is devoted to the derivation of our core result, 
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concerning PPT bound entangled states of minimal rank in 3 x 3 systems. As we 
mentioned above, they are all locally transformed projections onto the orthogonal 
complement of a subspace spanned by an orthogonal UPB. On the way to prove 
this result, we give numerous insights into how PPT states relate to product vectors 
in their kernels. That there is a strong relation between this subject and general 
Unextendible Product Bases becomes clear in Section 13.41 where we prove that 
there always is a minimal gUPB in the kernel of a non-separable PPT state of 
rank 4 in 3 x 3 systems, and we make much use of the general results of the first 
part of the paper. A common trait in all the following is also an idea of local 
SL (3, (D) (E) SL (3, (D) equivalence, explained with considerable care in Section [01 
The equivalence is natural both for gUPBs and positive-partial-transpose states, as 
well as in algebraic geometry when the Segre variety is concerned. It should also be 
noticed that we make a number of comments, related to the results of [14l[T^ and 
the question of atomicity of entanglement witnesses [T71IIH]- They are somewhat 
out of the main focus in this paper, but in our opinion they may prove valuable for 
the reader. 

Concerning the notation and language we use, only few things may need an 
explanation. Throughout the paper, {ci}^^^ denotes the canonical basis of (D". 
That is, Ci = [ ... 1 ... ] with 1 at the i-th position and zeros 
elsewhere. We denote by Keip, R (p) and r [p) the kernel, the range and the rank 
of a state p. The symbols Ti and T2 stand for the partial transpose with respect 
to the first and the second subsystem, respectively. For example, p^^ is a state 
p, partially transposed with respect to the second subsystem. We will be using 
the term "edge state" from [12] . A PPT edge state p is by definition a p which 
does not admit (p^ijj G R{p) such that (/)* (E)ip € R [p^^ ) . Equivalently, there is no 
<f><S)iJ e R{p) such that (f>^^/j* £ R (p^^). We say that a state p, acting on C"®C™, 
is supported on a fc x Z subspace when there exist a fc-dimensional subspace £1 of 
(D" and an /-dimensional subspace £2 of C" such that R (p) C £1 ® £2, while there 
exists no similar subspace C[ ® C'2 with dimC[ < I or dim £2 < k. We should 
also add that the notation E„_i^m-i, in accordance with is used for the Segre 
variety in p"™-i^ or less precisely, for the set of product vectors in (D" (81 (D™. We 
do not distinguish much between the two. 

2. General Unextendible Product Bases 
2.1. Basic results. Let us start with a definition of a general UPB. 

Definition 2.1. Take n,m S N. By a general Unextendible Product Basis, 

or a gUPB for short, we mean a set {0^ ® V'iliLi of product vectors in (D" ® C", 

< fc < mn, such that there is no product vector in span V'iliLi) ' 
orthogonal complement to the linear span of {(t)i '8)?/'j}i=i- 

In other words, a gUPB is a set of product vectors {(pi ® V'iliLi such that there 
is no product vector orthogonal to all of them. Yet another way of phrasing the 
definition is that the orthogonal complement to a gUPB is a Completely Entangled 
Subspace f20j, or CES for short. 

Remark 2.2. The definition of a gUPB can be trivially extended to a multipartite 
setting. 
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We know that gUPBs do exist. Any UPB consisting of orthogonal vectors is an 
example (concrete UPBs can be found e.g. in 2 ). We also know that for some 
spaces, no UPB consisting of orthogonal vectors can exist. For example, it has 
been noticed as early as in [1 that 2xn systems do not admit an orthogonal UPB, 
and a more general discussion of existence questions for orthogonal UPBs has been 
included in 8\ In the following, we show that gUPBs are much more common 
than the usual UPBs, and give a characterization of gUPBs of minimal number of 
elements. In Section 12.21 we discuss in more detail the 2 x n case, where much is 
known as a consequence of the recent work [16) . 

First, let us answer a question about the minimum number of elements in a 
gUPB in €" ® C™. 

Proposition 2.3. A set of vectors {(pi (E) V'iliLi ^ "E" (8) (C™ consisting of k < 
TO + 71 — 1 elements is not a generalized UPB. 

Proof. There exists a vector / e C" orthogonal to all the vectors (pi with i = 
1, . . . ,n — 1. Moreover, there exists a g G (D™ orthogonal to the vectors ipj with 
j — n, . . . ,k (because k — n < m). The vector f ^ g is orthogonal to all (pi (g) V'i for 
i^l,...,k. □ 

Proposition 2.4. A set of vectors {(pi (g) -0*}^t"~^ C (E) is a gllBP if and 
only if any n-tuple of vectors in {4>i}™i^ ^ consists of linearly independent vectors, 
as well as any m-tuple of vectors in {V'ilJlJi" ^• 

Proof. In order to prove necessity, assume that an n-tuple of vectors {0^, Yi^i is 
linearly dependent. Therefore there exists a vector / e (D" orthogonal to all of 
them. Vectors of the form f®g € (D"(8)(D'" with an arbitrary g g (D™ are orthogonal 
to all the vectors {p,^ (g) C C" ig) C™. Obviously, g can be chosen in such 
a way that / g is orthogonal to the remaining to — 1 elements of (D"* (because 
TO — 1 < TO = dim (D™). For a sufficiency proof, assume that f ^ g is orthogonal to 
(pi (S) ipi for i = l,...,TO + n— 1. The vector / can be orthogonal to at most n — 1 of 
(pi^s, whereas g cannot be orthogonal to more than to — 1 ipi's (remember the linear 
independence of n-tuples and TO-tuples, respectively). This gives a maximum of 
(n — 1) + (to — 1) vectors in {(pi (g) V'ilta"" orthogonal to / (g 5. Therefore / g) 5 
cannot be orthogonal to all the (pi (g ipi^s, the set {(pi g) V'ilS" ^ gUPB. □ 

As we already mentioned above, it is a trivial observation that the orthogonal 
complement to a gUPB is a Completely Entangled Subspace [3D]. However, it is 
more interesting to notice that the minimum number to + n — 1 of elements in 
a generalized UPB exactly matches a factor which appears in the formula D = 
mn — (to + n — 1) for the maximum dimension D of a CES in (D" g) [20U22j . 
This may seem a little bit puzzling at first, but can be easily explained. Consider 
a generic set of to + n — 2 vectors {(pi g) ipiY^^ ^ C C" g) (D™. It will have the 
property that any 71-tuple in {'/'i}™^""^ and any TO-tuple in {"i/'jIJI^"^^ ^re linearly 
independent. Therefore if there is a vector cf) (S) ip orthogonal to all the (pi g) ■0i's, 
we must have {cp, (pi) = for exactly n — 1 vectors (pi and {ip, tpj) = for exactly to 
TO — 1 vectors ipj. More precisely, there must exist AA, C {l,2,...,TO + n — 2} 
such that #M = n-l, #M = to-1, AfuM = {1, 2, . . . , to + n - 2} and {(p,(p^i) = 
OVigA^, as well as {tp^ipj) = OVjgjVi- Moreover, for any choice of an (n — l)-tuple 
Af of vectors (pi and a complementary (to — l)-tuple M. of vectors ipj, there exists 
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exactly one product vector cj) tp with = 0\/iej\f and {ip,ipj) = 0\fj(=M- 

Altogether, we have suitable pairs {Af,M.) of (n — 1) and (m — l)-tuples, 

which gives us exactly ("^^"7^) Pi'oduct vectors orthogonal to all the 0i (g) V'i's. To 
no surprise, this number equals the degree of the respective Segre variety I]„_i_m_i, 
cf. [m Lecture 18]. A key fact for us is that the number of product vectors is finite. 

If we assume that span ^{0^ (g) 7/;^}™^" and I]„_i.„i_i intersect generically 

enough [15l Example 18.5], we are lead to a conclusion that the projective dimension 
mn — 1 — (m + n — 2) of the orthogonal complement to the 0i ® V'i's must be bigger 
exactly by one than the projective dimension Z? — 1 of a maximal CES. Thus we get 
D = mn — (to + n — 2) — 1 = mn — (m + n — 1) , as expected. Note that this is only 
a rough explanation and not a rigorous proof. We also use it as a good opportunity 
to introduce and explain the formula for the degree of the Segre variety, which we 
are going implicitly refer to later. 

It is intuitively clear that the n-tuple (m-tuple) linear independence property of 
Proposition 12.41 is generic for random finite sets of vectors in (D" ((D™, resp.) and 
therefore minimal gUPBs should exist for any m, n (one can check it numerically, 
for example). However, it does not seem entirely trivial to prove this fact in a 
mathematically rigorous way. We do it in the following proposition. 

Proposition 2.5. For arbitrary n,m GlM, there exists a general UPB in (D" (K) C™ 
consisting of m + n — 1 vectors. 

Proof. According to Proposition 12.41 it is sufficient to prove that there exists a set 
of TO + 71— 1 vectors <j)i E (D" s.t. any n-tuple is linearly independent, and similarly 
a set of TO + n — 1 vectors ipi G C-'" with linearly independent TO-tuples. We prove 
a more general statement in the following lemma. 

Lemma 2.6. For any n € N, there exists an infinite sequence {ipi}^^ of vectors 
in (D" such that any n-tuple of them is linearly independent. 

Proof. Consider the family of vectors 

(1) 0. = [ 1 a, af . . . ] , 

where {cti}^^ is a sequence of numbers (real or complex) such that ai ^ aj for 
i j. Let us choose some {ii, 12, ... , in} G N s.t. ij ^ ik for j ^ k. The matrix 



(2) 



1 Qfj, af^ ... a^^ ^ 
1 a,, a\ ... a^-^ 



1 a,„ a|„ ... a--' 
is a Vandermonde matrix and it is well known that its determinant equals 
(3) W (a,^. -a^J. 

Therefore, the matrix is nonsingular for any choice of the n-tuple ^2, . . . , «„} G 
N. Equivalently, any n-tuple of vectors <j)i from equation ([T]) is linearly independent. 

□ 
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From Lemma 12.61 almost immediately arrive at the proposition. Choose 
m + n — 1 vectors (f)i as in equation ([T|) . In a similar fashion, take 

(4) = [ 1 A: ... /3r' ]: 



where again 7^ f3j for i ^ j. By Lemma 12.61 and Proposition 12. 4[ the set 

{(pi (g) V'ilti""^ constitutes a gUBP. □ 



It is natural to ask for a generalization of Proposition 12.41 for sets of product 
vectors consisting of more than m + n — 1 elements. We have the following 

Proposition 2.7. A set of vectors {4^ (g) C C" (g) C™ with N ^ m + n- 1 

is a glJBP if and only if for any A/", TW C N such that Af D M — and Af U M = 
{1,2, . . . ,N}, at least one of the sets of vectors {4'i}i£j\f '^i^d {V'iljgTvi spans the 
entire corresponding vector space ( (D" or (D™, resp.). 

Proof. Let us first prove necessity. Assume that the vectors {(f>i (g V'ili^i constitute 
a gUPB and choose some A/", C N as in the statement of the proposition. If 
neither of the sets {4'i\i£j\f {V'jIjgTvt spans the respective vector space, there 
exist / e (D" and g € C" such that = and (gjV'j) = for all i e 

J\f and j £ M. Because of the condition U M — {1, 2, . . . , iV}, we clearly 
have {f <E) g,(j>i (E) ipi) = for i — 1,2,...,N. This clearly contradicts the fact 
that the vectors g) tpi constitute a gUPB. In order to show sufficiency, assume 
that / g) g e C" g) C" is such that (/ g) 5, g) Vi) = for all i = 1, 2, ... , iV. 
Define A// := {i\{f,(f>i) = 0} and Mg := {j\{g,ipj} = 0}. Clearly, we must have 
TV/ U Mg = {1,2,..., N}. Thus it is possible to choose AA C A// and M C Mg 
such that AA U = {1, 2, . . . , iV} and A/" n = 0. By the very definition of A/> 
and Mg, we have (/, (pi) — for alH £ A/" and {g, tpj) for j G M. But according to 
the assumptions of the theorem, this is only possible if / and g are equal to zero. 
Thus {0j g) VJj}^! is a gUPB. □ 

Remark 2.8. For N = m + n — I, Proposition \2. 7| is the same as Proposition \2.4\ 
which should he expected. 

Remark 2.9. For a given set of vectors {pi g) N ^ m + n— 1, the condition 

of Proposition [^?7| can easily be checked, e.g. using a computer algebra program. 

Certain characterisations of gUPBs were earlier obtained in [lOj , but the above 
results were, rather surprisingly, never explicitly given. They can also easily be 
generalized to a multipartite setting. 

Proposition 2.10. A set of vectors {(p] g) 0f g) . . . g) 0-}^^ C C"ig)C""g). . . g)C"' 
with N ^ X]i=i ni + l~ \ is a gUBP if and only if for any Mi , . . . , A// C N such that 
Mi n Mj = for all i ^ i and Ij'^j^ Mi = {1,2,..., N{, at least one of the sets of 
vectors {(p) } j^j^j-. ; * = 1j 2, . . . , /, spans the entire corresponding vector space (D"* . 



Proof. Follows the same lines as the proof of Proposition 12.71 and will be omitted 
here. □ 



Note that a gUPB consisting of more than m + n— 1 elements does not have to 
contain a minimal one. This can be seen in the following example. 
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Example 2.11. Choose the following five vectors in (D^ ® C"^, 

(5) {ei (g) 61,62 (E) 62, (6i + 62) (E) 63,62 (E) (ei + 62 + 63) ,61 ® {ei + 2e2 - 63)} . 

It is a simple exercise to check that the set ([5|) satisfies the assumptions of 
Proposition l2.71 but neither of its four-element subsets is a minimal gUPB. However, 
it is still an open question whether there exists a minimal gUBP in the linear span 
of the vectors ([5]) . In fact it can easily be checked that the following four vectors 
belong to the linear span of ([S]) and they constitute a minimal gUPB. 

(6) (ei - 62) «) (61 - 262 -f 63) (61 -I- 62) ® (61 - 262 + 263) 



Thus, by giving Example 12.111 we have not proved whether there exist a linear 
subspace complementary to a Completely Entangled Subspace 20 and without a 
minimal gUPB in it or not. This seems to be an interesting question on its own. 

2.2. The simplest case. In the following, we are going to shortly discuss the 2 x n 
case, where much about generalized Unextendible Product Bases follows from the 
recent work [16] . Let us however begin with something which has been known for 
much longer time [T]. 

Orthogonal UPBs in 2 x n systems do not exist. In Ij, it followed from a local 
measurement protocol, but a direct proof can also be given, and it seems rather 
instructive to provide it here. 

Proposition 2.12. Let n € 'N be an arbitrary positive integer. A set of mutually 
orthogonal product vectors {(pi (E n C (D^ C" either spans the whole 

(D^ E) C", or it admits a product vector f E) g G (D^ (D", orthogonal to all the 

ig) V'j 's. 

Proof. For an inductive proof, let us first consider the case n = I. By Proposition 
12. 4[ a set of product vectors {(pi E) ei}-^-^ ^ in (D^ ® (D^ is a gUPB if and only if 
the vectors (pi span (D^. Equivalently, (pi E) ei's span (D^ E) (D^. 

Now assume that the proposition holds for all n = 1, 2, ... ,p — 1, p > 1 and 
consider a gUPB {(pi E) fpi}i^i ^ in (D^ (g) (D^, consisting of mutually orthogonal 
vectors. Define X :— {i\(f>i ~ (pi}, J :— {j\(pj^(p\\, and K, as the complement of 
lU J' in {1, 2, . . . , N{. By choosing M — Tva. Proposition 12. 71 we see that the vec- 
tors {V'jljgjuA: span (D^. By taking N = J , vfe arrive at a similar conclusion for 
{^iiiGiuK- Moreover, the vectors within the group {V'Oiei l^ave to be orthogonal, 
and similarly for {ipj}.^j. Both sets {^i}i^x ^"^"^ {'4'j}j^j also need to be orthogo- 
nal to {i/'/c}j,gx;- Taking all these facts about {V'iliei and {V'jljgj- together, we have 



span -g^:) = spa.n [{ipj}-^jj = span ({V'fel^.g^) , as weU as #1 = In 
the case #1 = #^7 = p, we have N — 2p, JC = 9, and the vectors {(pi E) fpi}i^i ^ 
span C^gjCP. If /C 7^ 0, {(pk E) V'felfeeK is a gUPB in C2(gspan {{ip k} kfzjc) , "^i^^ch is 
equivalent to E) C" for some n < p. By our inductive assumption, {(pk (E ^k}k£K 



(7) 



61 ® (61 - 662 + 863) 



(61 + 262) ® (61 - 462) 
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Let us now state the following lemma, which is simply an amalgam of Lemmas 1 
and 2 in a paper by other authors jl6j . 

Lemma 2.13. Let V be a Completely Entangled Subspace in (D^<8)(D". The product 
vectors (j) ®ip orthogonal to V span V-^ , while their partial conjugates (ff ® span 

We immediately get the following 

Corollary 2.14. Let y C (K) C" be a CES of dimension n - k, k ^ 1. The 
orthogonal complement admits a general Unextendible Product Basis, consisting 
of 2n — {n — k) = n + k elements. 

Proof. By Lemma [2. 131 we know that it is possible to find a basis {(f>i (g) of 
V^. Since is a CES, (g) '4'^}i^^ must be a gUPB. □ 

Let us also point out some algebraic geometry content of the question studied in 
Lemmas 1 and 2 of [16 . 

Remark 2.15. Let V be an {n — 1)- dimensional CES in (D^ (E) (D". As proved in 
Lemma 1 of I16f . the family of product vectors in has the form 

(8) (l,a)(»A(l,a,...,a"-i) 

with a g (D and A a non-singular n x n matrix. Such a family is a rational normal 
curve fl3\ Example 18.8], as a general intersection of the Segre variety with 
an n-plane in p2"-i js_ 

Proof. Obviously, vectors of the form ([5]) constitute a rational normal curve, cf. |15[ 
Example 18.8]. A general intersection of with an n-plane in p2"-i -^^111 be 

nondegenerate and transverse. By Bezout's Theorem [T5J Theorem 18.3], it will be 
a nondegenate curve of degree n in an n-plane. It is so because the degree of I]i^„_i 
equals n. Consequently, the curve must be a rational normal curve, cf. Proposition 
18.9 of [T5]. Thus we have explained our remark. However, it seems instructive to 
show why the intersection is transverse and nondegenerate, without repeating the 
whole argument of Lemma 1 in [16] . Let us first discuss the issue of transversity. 
Take a, (j)(E)tp in . As pointed out in fl6] , V is spanned by n — 1 product vectors 
^'i of Schmidt rank 2, with some additional conditions they must fulfill. They can 
always be written as <^ (g) i/'o + ® "01; where (0, (^^) — and (0^,0^) — 1, 
which uniquely defines the vector (j)-^ C (D^, as well as forces tpg and ip\ to be 

uniquely defined. Moreover, we know that {V'o}"— ^^'^ {^i}"— must be linearly 
independent. Otherwise, it is would be possible to find a product vector in V. As 
a conequence, tp is uniquely determined by the set of equations (^I^i, ® -0) = as 
the only vector in (D" that fulfills ('0°, '0) for alH = 1, 2, . . . , n — 1. Motivated by 
the description of the tangent space T^^^ (Si „_i) to the Segre variety at 4>® tl), 
provided in Lemma 13.51 from Section 13.31 let us consider vectors of the form (g) 
ip' + (j)' (g that belong to V-^. We can always assume 0' = 0-"-. After we do so, ^p' 
is determined by the set of equations (vPi, </> €5 0' -I- 0-*^ (X) "0) = Vi=i_...^„_i, which 
gives us (00,0') = — ('0i,'0): « = 1, 2, . . . , n — 1. As we mentioned above, the 
vectors {4'o}^^i i ^ ^^^^ {''/'iliLi^' linearly independent, and -0 is the only 
vector perpendicular to all the ipQ^s. Therefore the equations (0O' V'') = — (0i;0) 
have, up to a trivial scaling factor, a one-parameter family of solutions A'0+'^, where 
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^ is not proportional to ip and itself a solution to {ipQ, £,) — — (V'lj Therefore the 
intersection of T^®^ (I]i^„_i) and has dimension 2, or the projective dimension 
1. By a dimension counting argument practically identical to that presented in 
Lemma [3.61 we see that V'^ and I]i.„_i intersect transversely. In order to prove that 
the intersection is a nondegenerate curve, it is sufficient to show that it spans V-^. 
Obviously, this is one of the statements of the Lemma 12.131 but let us show it in a 
somewhat more abstract way. First of all, it is easy to show, by a little modification 
of the argument above, that the equations (^i, (X) -0), * = l,2,...,n — 1, give 
exactly one solution for ip when is fixed. Therefore we have a family of product 
vectors in V-^, (j) (E) ip {(j>). The precise form of the function ip (0), as in (j2.15p . is 
not important to us at the moment. Just as the authors of |16] did in the proof of 
their Lemma 1, it is fairly simple to show that the set {tp (0)}0g<g2 must span (D". 
Therefore it is possible to choose n vectors (pi G (D^ so that {ip {(pi)}^^^ is a basis 
of (D". Choose another cp, (pn+i, not equal to any of (pi with i = 1, . . . ,n. From the 
argument about transversity it is not difficult to see that ip {(pi) ^ {(pj) for (pi ^ (pj . 
A simple argument then shows that the family of product vectors {(pi ip {(pi)}"^i 
is linearly independent, and thus it must span V-^. This means nondegeneracy. In 
summary, the intersection of and I]i^„_i is transverse and gives a nondegenerate 
curve of degree n in an rt-plane, which must be a rational normal curve. This is in 
full compliance with the specific form ([S]) , obtained in . □ 

We should also warn the reader of a too simple conjecture, which may be ap- 
pealing because of the discussion above. 

Remark 2.16. It is not true that any CES admits a gUPB of its orthogonal com- 
plement. 

Proof. By [55] we know that Completely Entangled Subspaces are generic among 
all subspaces of fixed dimension k, as long as k is such that a fc-dimensional CES 
can exist at all. Let us look at the case of (D^ (8) (D^. The corresponding maximum 
dimension of a CES is 16 — 4 — 4+1 = 9, while the complementary dimension 
is 16 — 9 = 7. Both of these numbers are clearly smaller or equal 9. Choose a 
subspace V C (D^ ® (D^ of dimension 9 randomly. This is equivalent to choosing V-^ 
of dimension 7 randomly. Since both V and V'^ are chosen randomly, in the generic 
case neither of them contains a product vector, in accordance with the result of [23j . 
One can confirm this numerically. □ 

3. Non-separable positive-partial-transpose states of lowest rank 

The aim of the present section is to prove our main result, concerning positive- 
partial-transpose non-separable states of rank 4 in 3 x 3 systems. As indicated 
in [11], they all seem to be possible to locally transform to projections onto the 
orthogonal complement to a subspace spanned by an orthogonal UPB. Thus, there 
is strong numerical evidence that they are all locally equivalent to bound entangled 
states of the form discussed in [T]. In the following, we analytically prove that this 
is actually the case. Note that according to the results of [21], four is the minimal 
rank for an entangled PPT state. Therefore it is correct to say that our theorem 
concerns non-separable PPT states of lowest rank. Before we start with the proof, 
we need to discuss the concept of local equivalence, which will be fundamental to 
all that follows. 
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3.1. The concept of local equivalence. Numerous questions of physical or math- 
ematical origin need the proper identification of a group relevant to the problem 
in order to simplify the solution, or even to find it at all. The same is the case 
for the result we are going to obtain below. For PPT states, a natural group 
of symmetries should be of a product form, p i-^ (A (g) B)* p{A(E} B), because all 
such transformations preserve the property of being PPT. In physical terms, they 
preserve the splitting of a composite system into subsystems, which is a highly 
desirable property. The remaining question is, what group should A and B belong 
to. When the amount of entanglement between the two subsystems is in ques- 
tion, a natural choice is A and B in the Unitary or Special Unitary group. Such 
transformations cannot change any measure of entanglement. However, if the aim 
is to classify PPT states with respect to the property of being extreme, being an 
edge state [19] , or the number and dimensionalities spanned by the product vectors 
in their kernels or ranges, A and B should most naturally belong to the General 
Linear or Special Linear group. There is no essential difference between the two 
latter choices. Since we are not interested in positive scaling factors in front of 
the states, we choose to work with the Special Linear group. This was also the 
approach so successfully used by the authors of [TTJ[T3l[25]. We should remark 
that, while a PPT state is transformed according to p i— >■ (A (g) B)* p{A0 B), the 
product vectors in its kernel and its range undergo the following transformation, 
(j>^tp 1-^ {A^^ (g) B^^^ (t)®ip. Conversely, a transformation (/)(g)-0H> {A® B)(t)®'ip 
forces a change of p into p (^A~^ ® B~^^. It is these kind of trans- 
formations we will have in mind when we talk about "local equivalence", "local 
SL equivalence" or "SL (3, (D) SL (3, (D) equivalence" in the following sections. 
Any similar terms, even not listed here, will also refer to precisely the same situ- 
ation. Nevertheless, when product vectors in the kernel of a PPT state p are in 
question, it is more convenient to look at them as rays, points in the projective 
space. In such case, it is also more accurate to refer to the projectivisation of the 
group SL (3, (D) ® SL (3, (D), namely to PSL (3, (D) x PSL (3, C). In simple words, 
we may multiply vectors {4>i ® 'iPit4'2 ® ipiT ■ ■} C Kerp by arbitrary individual 
factors, and they will remain elements of the kernel of p. Together with the pre- 
viously introduced product SL (3, C) ® SL (3, (D) symmetries, we have a group of 
transformations that is most properly described as PSL (3, (D) x PSL (3, (D). Note 
that the use of this term is motivated mainly by the possibility to avoid excessive 
comments about constant factors in front of the product vectors in Kerp. We are 
fully legitimate to use the previously introduced name "local equivalence" also for 
the PSL (3, (D) x PSL (3, (D) transformations we just described because constant 
factors are completely irrelevant to p itself. 

The ultimate reason for using equivalences of the form described above will be 
the simplicity of our main result, a characterization theorem that we are going to 
obtain in Section 13.71 The equivalence classes under SL (3, (D) ® SL (3, (D) of non- 
separable PPT states of rank 4 in 3 x 3 systems turn out to be parametrized by 
just four real, positive numbers. Moreover, each class has a representative which 
is a projection onto a Completely Entangled Subspace complementary to a 3 x 3 
orthogonal UBP. This is a surprising result for which strong numerical evidence 
was provided in I 111 and later supported by certain analytical results of jl3j . 

3.2. Outline of the proof. The proof is not excessively complicated, but it needs 
a considerable amount of work. It also consists of a number of steps which do not 
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seem easy to merge. In order to simplify the reading, we start with a hst of building 
blocks. We will elaborate on each of them in the following sections. 

(1) The kernel of a rank four PPT state p must intersect the Segre variety in 
a transverse way. In particular, according to the Bezout's Theorem, the 
intersection must consist of exactly six points. 

(2) The product vectors in the kernel of a rank 4 PPT state in the 3x3 case 
span the kernel. As a result, they must be a generalized UPB. There cannot 
exist a product vector orthogonal to all of them. 

(3) A generalized UPB in the 3x3 case is locally equivalent to an orthogonal 
one if and only if certain invariants Si, . . . , S4, introduced in [11] . are all 
positive, possibly after the vectors are permuted. 

(4) A generalized UPB in a 3 x 3 system is contained in a kernel of some rank 
four PPT state if and only if the corresponding values of si,...,S4 are 
positive, possibly after the vectors are permuted. Moreover, in such case 
the PPT state in question is uniquely determined. 

The final conclusion from the facts mentioned in items (1) — (4) is that the only non- 
separable PPT states of rank 4 in 3 x 3 systems are local transforms of projections 
onto orthogonal complements of orthogonal pentagram-type Unextendible Product 
Bases. 

3.3. Product vectors in the kernel of a PPT state. The present section elab- 
orates on item (1) in the list given above and on related topics. Let us start with 
an elementary fact 

Lemma 3.1. A product vector is in the kernel of a PPT state p if and only 

if the partially conjugated states cf)* ® 'ip and (j) ® ijj* are in the kernels of p"^^ and 
p^'^ , respectively. 

Proof. It simply follows from the identity {(j) ^ ip , p {<j) = {(j)^ ip* ,p'^^ [<f) ® ip*)) 

and {4> ® , p {4> ® ip)) = {(j>* (S) ip , p'^^ {(j>* ^ip)), by the positivity of p, p^^ and 
p^\ □ 

In the above lemma, we did not assume anything about the dimensionality of 
the system. Neither we do it in the following. 

Lemma 3.2. Assume that a product vector cf) ® ^ is in the kernel of a PPT state 
p. In such case 

(9) {(p' (^^,p{4)®i!')) ^{(p®^',p{4)' =Q Vc^',^' 

Proof. Since p{(j) <^ip) — 0, we know from Lemma [3T] that p'^^ {(p* — 0, which 
obviously implies (</>'* (g) tp', p^^ {(p* ® ip)) ^ {(p<^ tp' , p {(p' ® -0)) = 0- This is the 
first equality in (|9|). The second one can be obtained in a similar way. □ 

Our next lemma applies specifically to the 3 x n case. For reference concerning 
edge PPT states, check l9|. 

Lemma 3.3. Assume that both (p® %p and cp' (g) tp, with (p, cp' in (C^ and tp in (D", 
(p (p' , belong to the kernel of a PPT state p, acting on ® (D'^. The state p 
is either supported on a 3 x [n ~ 1) or smaller subspace, or it can he written as 
p' + X\(t)" (E)0 for some A > 0, ^ € C", cp" £ linearly independent 

of (p and (p' , and a PPT state p' , supported on a 3 x {n — 1) or smaller subspace. 
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Moreover r (p') = r (p) — 1 and r 



{p')'^^j — r (p"^^) — 1- In a situation when the 



reduction is possible, the state p is not an edge PPT state. In particular, p is not 
an extreme and non-separable PPT state. 

Proof. Let us assume that the product states (j)' ® i}} and (/)2 (X" "0 belong to the 
kernel of p. Let A be an SL (3, (D) transformation that brings ei,e2 C (D^ to 0i 
and 02- A little inspection shows that Lemmas 1 and 2 of [24] can be applied 
to p :— [A®!) p{A® 1). Consequently, we see that either p is supported on 
a 3 X (n — 1) or smaller space, or the assertion of Lemma 2 of [53] tells us that 
p — pi + A |e3 (g) ^) (ea (g) ^1 for some ^ g (D", and moreover, pi is a PPT state 

supported on a 3 x (n — 1) or smaller subspace, with r (pi ) — r (p) — 1 and r ( j — 



r(p^i)-l. Wehavep= (A^i «) l)* p {A-'^ (g) l) = p' + A ^) ((/)"' ® where 
0"' = A^^es and p' = (A^^ (g) l)* pi (A^^ (g) l). The states p' and (p')^' stih have 
their ranks reduced by one with respect to the ranks of p and p'^^ , respectively. 
The subspaces on which they are supported are of the same type as for pi , hence 
3 X (n — 1) or smaller. The statement that p is not an edge state simply follows 
because \(j>"' g) ^) {(/)'" (E) £, \ is in i? (p) while its partial conjugation is in R (p'^^ ) . □ 

The following result reduces a more general case to the situation considered 
above. However, this time we assume n = 3. 

Lemma 3.4. Let cj) (E) ip G (C^ (E) (D^ be an element of a PPT state p, acting on 
(D'^g) (D'^. There cannot exist a nonzero vector ip'^ip' + (f>'(E4', with (j)' ^ (f> or ijj' ^ ip, 
in the kernel of p, unless one of the following is true: i) p = p' + X \ C ® £,) {C E) £,\ 
for A > 0, ^, C e (C^ and p' a PPT state supported on a 2 x '3 or smaller subspace 
with r (p') = r (p) — 1 and r (^(p')^^^ = r {(p)^^^ — 1 or ii) p is supported on a 2 x 3 
or smaller subspace itself. 

Proof. Assume that there is a state of the form ® -0' + 0' ® -0 in the kernel of p. 
This is equivalent to saying that {</> E) tp' + <j>' <E) t/j, p {(p E) ip' + 4>' E) = 0. The 
inner product factorizes as 

(0 «) 0', p (0 «) tA')) + (0 «) Ip', P (0' (E ll>)) + (0' «) -0, P (0 «) -0')) + (0' «> -0, P (0' ® -4^)) 

The two factors in the middle vanish according to Lemma [3.21 while the two remain- 
ing factors are nonnegative as a consequence of positivity of p. Therefore, the only 
possibility for the above expression to vanish is when {4> E) ip' , p {(/) E) ip')) = and 
(0' E>'4', p (0' (8) ■0)) = 0. This in turn means that p{(j)' E)ip) =0 and p{(j)E) ip') = 0. 
According to our assumptions, at least one of these equalities is nontrivial (i.e. 
0' 7^ or ijj' 7^ 0). Lemma [3731 can be applied. □ 

The importance of Lemma 13.41 is evident if we realize that the tangent space to 
the Segre variety, or to the set of product states at a point 4>E)ip, consists precisely 
of the vectors of the form considered above. We have 

Lemma 3.5. Elements of the tangent space to the Segre variety, or to the set of 
product vectors at a point (piE) ip, are of the form 




(10) 



(j) E) Ip' + (p' E) Ip 



with Ip' and 0' arbitrary. 
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Proof. A heuristic proof may consist in writing {(p + S<j>) ^{tp + Sip) « (p®ijj + S<p(S) 
ip + cpiSiSip, where the approximate equahty holds to the first order. A shghtly more 
rigorous one can be found in Example 14.16 of [TS]. □ 

Next, we specify the rank of p to be 4 and keep the assumption that p acts on 
(D'^ (8) (D^. Thus the kernel of p is of dimension 5, which is the smallest number d 
such that a d-dimensional linear subspace must intersect the set of product vectors 
in (D'^ (8) (D'^, cf. e.g. [50]. Following Lemmas 13.41 and 13.51 we can show that the 
nonempty intersection is generic in the sense of Bezout's theorem [TS| Theorem 
18.3] and thus it consists of exactly six points. 

Lemma 3.6. Let p be a non-separable PPT state of rank A acting on (D'^OC"^. The 
intersection between the respective Segre variety and the five- dimensional kernel of 
p is transverse at every point. There are exactly six product vectors in the kernel 
of P- 

Proof. The preceeding lemmas and the arguments included in their proofs imply 
directly. Let us take (p ^ ip G Kerp. As we mentioned above, such a vector exists 
[20l[22] by a dimensionality argument for projective varieties. We easily see from 
Lemma lSTSl that the dimension of the tangent space T^^^ (^2,2) to the Segre variety 
at 0(g)?/' is 5, and thus the projective dimension is 4. Being more explicit, any vector 
of the form (j)®\p' -\-(p' ®\p can be written as \4i®ip + Y^1^i £,iP®4'i+Y^j=i Cj<Pj®i^i 
where {</>, (/>!, (j)2\ and {ip, ipi, "02} are two sets of three linearly independent vectors 
in (D^ and Xi, (j are arbitrary complex coefficients. From Lemmas 13.41 and 13.51 we 
know that the only vector in the intersection of Kerp and T^,^^ (^2,2) is <p<E)ip itself. 
It must be so, because otherwise we could reduce the rank of p by subtracting a 
projection onto a product state. After the reduction, we would be left with a PPT 
state of rank 3. However, all such PPT states are separable according to |24j . 
and p would have to be separable as well. The other option is that p could be 
supported on (D^ (E) C'^ or even a less dimensional space itself. But then it is well- 
known that p is separable as as consequence of being PPT [35]. In either case, 
we get a contradiction with the assumption that p is non-separable. Therefore, 
p (S) ip must be, up to a scalar factor, the only element of the intersection between 
Kerp and T^^^ (^2,2)- Consequently, the dimension of Kerp-|- T^^^ (^2,2) equals 
5-1-5—1 = 9, and the projective dimension is 9 — 1 = 8. This equals the projective 
dimension of T^j^^ (-^P*)' simpler, the dimension of the projective space P®. In 
other words, Kerp and T^^^ {^2,2) span T^ig,^ i^^)' '^tiich is equivalent to saying 
that the intersection between Kerp and the Segre variety is transverse at (p iSi ip. 
Since we did not make any additional assumptions about p ® ip apart from that 
it belongs to the intersection, we see that the intersection is transverse at every 
point. Therefore Bezout's theorem applies. The fact that there are exactly six 
points in the intersection follows because the degree of the Segre variety 112,2 is 
sbc [151 Example 18.15]. ' □ 

In summary, in the present section we have shown that a non-separable rank 4 
PPT state in a 3 X 3 system must have exactly six vectors in its kernel. This is in 
full argeement with an assertion of [14]. It should be noticed that, as a part of the 
proof of the above lemma, we have shown that non-separable PPT states of rank 4 
in 3 X 3 systems are edge states. Thus, Lemmas 13.31 andl 3. 41 can be directly applied. 
We will frequently use them in the following section. 
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3.4. Product vectors in the kernel must be a gUPB. We already know that 
the number of product vectors in the kernel of a rank 4 non-separable PPT state of 
a 3 X 3 system is six. In the following, we discuss more specific properties of the set 
of six product vectors. Let us denote them with 0^ ® -0^, i = 1, 2, . . . , 6. It turns out 
that, up to local equivalence, five of them can always be brought to a special form, 
which has only four real parameters, the numbers Si, . . . , S4 introduced in [11'. It 
then follows that the vectors (pi if they belong to the kernel of a rank 4 PPT 

state, must span a five-dimensional subspace. Thus they span the kernel. 

In order to prove our assertion, first observe that (pi ^ 4>j for i j, and thus 
they must span at least a two-dimensional subspace of C"^. Similarly for the ?A's. 
Let us assume first that one of the set {(f>i}^i^i or spans a two-dimensional 

subspace. We may, for example, try to assume this about {<pi}^^i- Up to PSL (3, (D) 
transformations, we have 

"101111 

(11) [ (pi (f>2 (f>3 (t>4 4'5 (pe ] = I I p q r 

L 

where p, q, r are all different and different from and 1. When writing ()11|) . we 
used the fact that there is no pair of identical vectors in {4>i}i=i- Up to local 
transformations, we have ■01 = ei and "02 = 62. As for the other vectors ip, we use 
the following notation, ^i = [ ipu 1^21 tpsi ], i = 3,4,5,6. We also introduce 
coordinates w'-' for general vectors w — ^ ^3 <D^ €5 (D"^. Our aim is 

to show that there exists a linear combination of the vectors (pi ipi of the form 
(j> ^ tp' + (j)' ^ ip from Lemma This will lead us to a contradiction and show that 
0i's cannot be as in (ITT|) . and must span (D'^. An analogous conclusion for 0's will 
be immediate. 

Let us first observe that ^p3i 7^ for all i e {3,4,5,6}. Otherwise, we would 
have three product vectors supported on span (61,62) (8) span (61, 62). Up to local 
equivalence, they would be of the form 61 (X) 61, 62 (81 62 and (61 -I- 62) (8 (61 -I- 62). 
In such case, 61 62 -|- 62 8) 61 = (61 -I- 62) (8 (61 + 62) — 61 8) 61 — 62 (8 62 would 
be in the kernel of p, which contradicts Lemma 13.41 Therefore we must have 
031 7^ for all i. Let us choose a and /3 so that aV'33 + Ppi'sA — 0. The vector 
a03 8) -03 + P4>A ® 04 has a vanishing coordinate = 0-033 + l^pipsA and a non- 
vanishing coordinate uj^^ = aip^s + I3ip34 (remember that p 7^ 1). By subtracting 
62 €5 62 times aip23 + Pq'ipai, we can cancel the coordinate, and similarly cancel 
by subtracting aipis + jSipu times 61 (g) ei. In the end, we see that a vector of 
the form cj^^62 ® ei + ui^^ei ^ 62 + ^-'^^61 (8 63 with uj-^^ 7^ is in the kernel of p. 
But this contradicts Lemma 13.41 In summary, the vectors 0^ cannot be brought to 
the form (fTT|) . or in other words, they span (D^. Obviously, the same is true for the 
set {V'ili^i- A more careful analysis of the above argument leads to even stronger 
conclusions. Firstly, an assumption that there exist three vectors 0i €5 ipi supported 
on a 2 X 2 dimensional subspace lead us to a contradiction. Therefore we have the 
following 

Lemma 3.7. Let {(pi (g) 0'i}^=i defined as above. For any triple |0.i^. (g) ipi. C 

{(pi (g) '0i}i=iJ least one of the sets of vectors {(pi^ •^P'^^'S '^^ ■ 

Moreover, we only needed four product vectors with 0's as in pT|) to arrive at a 
contradiction with Lemma 13.41 As a consequence, we have 
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Lemma 3.8. For any quadruple {^^^ ® ''l^ij") j-i ^ both the sets of 

vectors {<t>ij^'^j^i o^^d {ipij}*^-^ span (D^ . 



As an immediate consequence of Lemma 13.71 there exists a set of three Un- 
early independent vectors in {4>i}^^i- With no loss of generahty, we may assume 
that {01,02,06} is a hnearly independent set. After a PSL (3, (D) transformation, 
01 = ei, 02 = 62 and (f>6 = S3- There are in principle two possibilities concerning 
the remaining vectors 03 , 04 and 05 . Either one of them is of the form [ x y z ] 
with xyz 7^ 0, or all of them have exactly one coordinate equal to zero. Two vanish- 
ing coordinates in a single vector cannot occur because there is no pair of identical 
vectors among 0i, . . . , 06. Moreover, according to Lemma 13. 8( the zeros must oc- 
cur in different places in 03, 04 and 05. Up to PSL (3, (D) transformations and 
permuting the vectors, we may assume that 03 = [2: l],04=[O 1 z], 
05 = [ 1 y ] with X, y, z all different from 0. But then, write the coordinate 
matrix for {01,02,03,04}, 



(12) 



10x0 
10 1 
1 z 



It is easy to check that all the 3x3 minors in (|12p are non-vanishing. In other 
words, any triple of vectors in {0i, 02, 03, 04} spans (C^. The corresponding vectors 
"01, tp2, V'3 s-nd '04 may or may not have all triples linearly independent. If all 
the triples span (C^, we can simultaneously, by using a PSL(3,(D) x PSL (3, (D) 
transformation, bring {01,02,03,04} and {ipi,ip2,'4'3,ipi} to the form 



(13) 



[ 01 02 03 04 ] = [ tpl i>2 i>3 1p4 



10 1 
10 1 
11 



By adding a fifth product vector, say 05 (S) 1^5, we get, up to local transformation 
and relabelling the vectors 0^ (g) -0^, 



(14) 



-01 ^2 V'3 ^04 V'5 



10 11 
1 1 p 
1 1 g 



10 11 
1 1 r 
11s 



where p, q, r, s are some complex numbers. We should remark that the possibility 
to have 1 in the first coordinate of 05 and follows because there must exist i S 
{1,2,3} such that 0i5'0i5 7^ 0. Otherwise, 05 or -05 would have to be proportional 
to Ci for some i e {1, 2, 3}. 

If not all triples in {^pl,^p2^tp3^''P'l} s-re linearly independent, it is still possible, 
according to Lemma l5^ to find a linearly independent triple among them. Without 
loss of generality, we may assume that the triple is {?0i, -02, "03}- By an identical 
argument as for the 0's, we know that there is a vector -0^, i € {5,6} such that 
{■01, '02, V'Sj ''pi} have all triples linearly independent. Without loss of generality, we 
may assume that 0i = 05. This time, a local transformation and possible relabelling 
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brings the product vectors 0i (8) V'i with i — 1,2,. ..,5 to the form 

10 11" 
1 1 p 
1 1 g 

10 11 ■ 

1 r 1 
1 s 1 _ 

To make a final touch to this section, we need to show that product vectors of 
the form (jl4p or (jl5p are hnearly independent if they are five in number, and thus 
they span the five-dimensional kernel of p. We will also show that they constitute a 
minimal gUPB, and that the parameters p, q, r, s have to be real when the vectors 
are in the kernel of a PPT state. 

Let us use the notation [ ] for 

vectors w = j oj^^ei (8) Cj in (D'^ (g) C"^. In the case p4|). we have 



(16) 



" 01 a 






" 1 


























02 


5 ^^2 
















1 














03 a 


5 V'3 




























1 


04 


5 ^4 




1 


1 


1 


1 


1 


1 


1 


1 


1 


. 05 


5 '05 . 




1 


r 


s 


p 


pr 


ps 


q 


qr 


qs 



In the case ([T5]) , the coordinates of the product vectors are the following. 



(17) 



" 01 <i 


5 -01 




" 1 


























02 


5 V'2 
















1 














03 « 


3 V'3 




























1 


04 ^ 


5 V'4 




1 


r 


s 


1 


r 


s 


1 


r 


s 


. 05 <i 


?'V'5 . 




1 


1 


1 


p 


P 


P 


q 


q 


q 



It is an elementary exercise to check that the matrices on the right-hand side of 
(|16p and (|17p are of rank 5 for all choices of p,q,r,s, with the only exception 
ofp = q = r = s = l. But the last possibility is excluded because it implies 

04 (g V'4 = 05 (X) V'5- 

Next, we can show that the vectors 0^ tpi with i = 1,2, ... ,5, chosen as above, 
constitute a general Unextendible Product Basis. In order to prove it, let us first 
show that the rank of p^^ has to be 4. 

Proposition 3.9. Let p be a non-separable PPT state of rank 4 acting on <D^(E)(D^ . 
The rank of p^^ is also 4. 

Proof. If p is non-separable, we know by the above argument that the product 
vectors {0^ (S) V'iliLi the kernel of p span a five-dimensional subspace, which is 
the kernel itself. Moreover, five of them are, up to local transformations, of the form 
([T^ or But this implies that the corresponding product vectors in the kernel 
of , which are 0* (g) tpi according to Lemma 13.11 can also be brought to the form 
(HH) or P^S)) . To be more explicit, if a local transformation A^B brings the vectors 
0i (g) tpi with i = 1, 2, . . . , 5 to the form or (fT5|) . A* i? does the same to the 
partial conjugations cj)* The only difference is that p and q change into p* and 

q* in ()14p or ()15p . But this does not change the conclusion about the dimensionality 
of the subspace spanned by vectors of the form (fT4| or (fTS]) . As a consequence, 
the product vectors in the kernel of p'^^ span at least a five-dimensional subspace. 



(15) 



V'l V'2 V'3 V'4 V'5 
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Thus the kernel of is at least five-dimensional. If it had higher dimension, the 
rank of p^^ would be lower or equal 3, which is, according to [24], impossible for 
non-separable p. Therefore, the dimension of the kernel equals 5, and the rank of 



There exist separable states p of rank 4 in 3 x 3 systems that have the rank of 
p"^^ different from 4. However, our next proposition shows that if p is supported 
on (g) and it cannot be written as p' -|- A |C (8) (C ® ^1 with A > and p' 
supported on a 2 x 2 subspace (cf. Figure 4 in [T31), the rank of p^^ is also 4. 

Proposition 3.10. Let p he a separable state of rank 4 supported on (E) (D^ , 
which cannot be written as p' + \ \C ® C) (C "8) ^\ with A > and p' supported on a 
2x2 subspace of (D^ (g) (D^. The rank of p^^ is also 4. 

Proof. First, we should remark that r {p'^^) = r {p^^)- This fact will be important 
for some parts of the proof, although never explicitly refered to. The main idea 
that we are going to use is that the argument preceeding formulas ([T^ and ([15]) 
works for separable states as well, provided that they cannot be reduced according 
to Lemma 1331 In other words, the argument works when the kernel of a PPT state 
in question does intersect the Segre variety in a transverse way. Thus, if a reduction 
according to Lemma 13.41 is not possible for a separable state p, we have vectors of 
the form (jl4p or (|15p in Kerp, and they span a five-dimensional space. This is also 
the dimensionality of the subspace spanned by their partial conjugates, which are 
in Kerp^i . Therefore the rank of p^^ is not bigger than 4. If it was less than four, 
the intersection between Kerp-^^ and the Segre variety ^2.2 would be more than 
zcro-dimensional, according to the Projective Dimension Theorem [27, Theorem 
7.2]. But this contradicts the fact that there are only a finite number of product 
vectors in Kerp-^^ (equal to (p* ^tp for all (p^tp £ Kerp). Therefore the rank of p^^ 
has to be 4 when Kerp intersects the Segre variety transversely. If not, we know 
from Lemmas and 13.41 and 13.51 that there are two options: i) it is possible to write 
p as p' + X\C, iC ^1) where A and p' is a rank 3 PPT state supported on a 



2 X 3 or smaller subspace of (D^ C^, with r (p') = 3 and r ( (p') M = d - 1, ii) p 



is supported on a 2 x 3 or smaller subspace itself. Option ii) is excluded because 
of the assumption of p supported on (D'^ g) (D"^. Our aim in the following will be to 



show that r (p') = r I (p') ^ J unless p' is supported on a 2 x 2 subspace, which is 



precisely the second possibility we allow in the proposition. First, observe that if p' 
is supported on a 2 x 3 subspace, we can use an analogue of Lemma l3.4l Either we 
have p' = p" + A' ](' (g) (C' <8) C'l where A' > and p" is supported on a 2 x 2, 1 x 3 



or 1 X 2 subspace, r (p") = 2 and r ( (p") ^) — t [ (p') M — 1, or Kerp' intersects the 



respective Segre variety Ei_2 transversely. In the latter case, by Bezout's Theorem 
the 3-dimensional kernel of p' has precisely three product vectors in it. Actually, 
we can repeat the argument preceeding Lemmas 13.71 and 13.81 to conclude that the 
product vectors in Kerp' have to be locally equivalent to 




□ 






01 (f>2 03 



1 1 
1 1 



(18) 



Ipl 1p2 1p3 



1 
1 
1 
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Obviously, these vectors span the kernel. We see that there are, within the 2 x 



3 subspace, only three product vectors in R{p') 
equivalent to 



(Kerp') . They are locally 



(19) 



Cl C2 Cs 



6 6 6 



1 

1 



1 



-1 



1 
1 
1 



Since p' is separable and of rank 3, it must be locally equivalent to a convex sum of 
projections onto the vectors C,i ® in , which implies that p^^ is an analogous 
sum of projections onto Q ® But Q ® = Q £,i if the product vectors 

are as in p9)) . Therefore R{p') = r i{p')^^ \ , which implies R{p) = r(p-^i), as 



expected. This proves our assertion for p' supported on a 2 x 3 subspace with Kerp' 
that intersects the corresponding Segre variety Si_2 transversely. For the other 
nontrivial cases, we can have p" separable and of rank 2, supported on a 2 x 2 
subspace. There is also the trivial case of p" supported onalx2orlx3 subpace, 
in which the equality r [p") = r ^(p")"^^) clearly holds, and it implies equality of 
ranks of p and . 

In the case of p" supported on a 2 x 2 subspace, we can repeat the argument 
with transverse intersections. Either p" can be reduced once again, in which case 
it turns out to be equal to A'" |C"' (C" ® + A" |C" ® C) (C" C'\ with 
A" > 0, A'" > and ® C not proportional to (" ® C' , or Kerp" must intersect 
the respective Segre variety Si^i in a transverse way. The first possibility clearly 

gives us T {p") = 2 — r (^{p")^^^- The latter implies, by Bezout's Theorem, that 

there are exactly two product vectors in Kerp". Similarly as for (|18|) . we can prove 
that the two product vectors must be locally equivalent to ei ei and 62 (8) £2. 
Clearly, they span the kernel of p" and there are only two product vectors, locally 
equivalent to ei (g) 62 and 62 (81 ei, in R{p"). But p" is separable and of rank 2. 
Therefore it must be locally equivalent to a convex sum of projections onto these 
two vectors. Accordingly, (p")'^^ is locally equivalent to a sum of two projections 
onto product vectors, which are e* 62 and 62 (8 ei, actually equal to ei (8) 62 and 

628)61. This clearly implies r (p") = r (j^p")^^^ and the equality between the ranks 

of p and p^^ clearly follows. □ 



Remark 3.11. The two propositions above explain why PPT states of ranks (4, n), 
n A should not be expected to appear in the upper part of Table II in \14^ . They 
do exist, but they are always separable and of a rather special form. 

It is useful to formulate the following 

Corollary 3.12. All rank 4 non-separable PPT states p in 3 x 3 systems are edge 
states. 



Proof. If some non-separable p of rank 4 had a product vector (j) <^ ^ in its range, 
and the partial conjugated vector cj)* ^tp was in the range of p^^ , we could diminish 
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the rank of p or p^^ by subtracting X\(t) ® {(}) ^ where 

(20) A = min|(0® (0® V-, ('^® ^)) ^} ' 

cf. [19]. In such case, p could be written as p = p' + A |(/) (g) i/") V"! with p' PPT 
and of rank 3 or with p^^ of rank 3. But this imphes, by [24], that p' would have to 
be separable. This further implies separability of p, which is a contradiction. □ 

At this point, we can easily prove that the vectors <f>i ® ipi in the kernel of 
a non-separable p of rank 4, chosen as in (fH)) or p3)) . constitute a generalized 
Unextendible Product Basis. If there was a product vector <^ (g) -0 orthogonal to all 
of them, it would be an element of the range of p. From the proof of Proposition 
13.91 we know that the partially conjugated vectors cj)* (E) ipi span the kernel of p^^ . 
Since {(p i^i) = = (0* (g) -0, (^* (g) ipi) for all i, we see that (g) V is in the 

range of p^^ , (Kerp-^^ ) . Therefore we have a product vector (g -0 in the range 
of p such that its partial conjugation is in the range p^^ . In other words, p is not 
an edge state. But this contradicts Corollarv l3.12l and therefore cannot happen. In 
this way, we have proved the following. 

Proposition 3.13. Let p be a rank 4 non-separable PPT state m a 3 x 3 system. 
The six vectors in the kernel of p constitute a generalized UPB. There is a subset 
of five of them that constitutes a minimal gUPB in the sense of Proposition \^^\ 

Proof. Most of the proof has already been explained above. We only need to com- 
ment on the fact that five of the product vectors constitute a minimal gUPB. It 
must be so because the five vectors we brought to the form (fT4|) or (fT5| span the 
kernel of p, and the orthogonal complement to the kernel has no product vector 
in it. Thus, the five vectors are a gUPB of Kerp, which is minimal according to 
Proposition l2.4[ because m + n— 1 = 5 for m = rt = 3. □ 

By Proposition l2 .41 we know that a minimal gUPB {(pi g) V'i}^=i ^^s the property 
that all triples in {4>i}^i^i and in ^i'^ linearly independent. In such case, 

the forms (fH)) and ([T5|) are locally equivalent, and we may choose to work with 
only one of them. In the sequel, we prefer to assume the form (|14p of the product 
vectors, which is in agreement with the convention used in [T?| . Our next step is 
to prove that the parameters p, q, r and s in (|14p must be real if the corresponding 
product vectors belong to the kernel of a rank 4 PPT state in the 3x3 case. This 
is not of much use here, but will prove to be important in Section [3.61 

We know from Lemma [3 . 61 that there are exactly six product vectors in the kernel 
of p, while we have only five of them in (jl4[) . and we know that they span the kernel. 
Consequently, the sixth vector is a linear combination of the other five ones, 

5 

(21) 0g g) -i/ig = ^ Ai^i g) Ipi 

i=l 

Note that explicit formulas for the sixth vector can be found in |131 Section 5.2]. 
Interestingly, since (pe0ipe S Kerp, we know from Lemma [X^ that (pl^^pe G Kerp 
is in the kernel of p^^. However, the vectors cp* (g ipi with i — 1, 2, . . . , 5 are also 
there and moreover, since they are, up to local equivalence, of the form (|14p with 
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p and q complex conjugated, we already know that they span Kerp^ . Thus 

5 

(22) ^*e(E)i,6=Y.^^^* <»^^, 

4=1 

where the coefficients are in principle not related to the A^'s in (pij) . However, 
we can already see at this point that it may be very difficult to simultaneously 
satisfy equations (|2T|) and ((22|) . if we do not assume that (f)i = (j>* for all i. In the 
latter case, one can obviously choose — A^. Our aim in the following will be to 
show that ~ Xi is the only possible choice. By projecting (PT|) onto the first, the 



second and the third coordinate in the first subsystem, we get 

(23) Ai?Ai + Xi-ipi + AsV's = (1)16 i>6 

(24) X2i'2 + XilpA + pXbi^b = (t>26lp6 

(25) AsV's + >^4tp4 + = </'36 V'6 

where {'/'ieji^i are coordinates of (pe- Similarly, from (j22p we get 

(26) 6^1 + C4V'4 + ^5^5 = 0^6 V'6 

(27) 6V'2+C4V'4+P*C5l/'5 = 026 V'6 

(28) ^V-a + C4V'4 + 9*^5^5 = 036 V'6 



Let us note that the triples {t/ji, V'4, ^5}, {^'2,^4,^5}, {V'37^47V'5} a-U consist of 
linearly independent vectors, according to Proposition 13 . 13l This implies that each 
of the formulas (|23p - (|28|) gives exactly one solution for the coefficients Xi or 
which it contains. For one of the consequences, all the coefficients ■0^6 must be 
non-vanishing. Two of them cannot vanish, because 06 proportional to any of 0i 
with i = 1, 2, 3 would contradict 06 V'e 7^ 04 ® tpi or Lemma [3.31 To see this, let 
us assume that one of them vanishes, e.g. 036 — 0. In such case, equation (j25p 
implies A3 = A4 = A5 = 0, where we used the fact that <? ^ 0. Hence ([23]) and 
((24)) reduce to 0i6V'6 = Ai-0i and 026'06 — A2'02- But neither of these equalities 
can hold, since 0ig 7^ and 026 7^ and V'e proportional to V'l or -02 contradicts 
Lemma 13.31 Thus our assumption 036 = must have been false. By repeating the 
same argument for 0i6 and 026, we arrive at 016 026 036 7^ 0- Let us also notice that 
necessarily A4 7^ and ^4 ^ 0. We cannot have, for example ^iV'i + CsV's = 016^06 
and ^2V'2 +P*£,5'<p5 = 026 V'6 since the only vector in the intersection of span {ipi, ip^) 
and span (-02 7 "05) is V'Si and we know that -06 7^ 0'5 by Lemma 13.31 or simply by 
06 V'6 7^ 04 "04- In a similar way, one obtains A5 7^ and ^5 7^ 0. With such 
amount of knowledge, we can easily prove the expected result. 

Proposition 3.14. Let 0^ (g) ipi for i = 1, 2, . . . , 5 be product vectors of the form 
(|14p in the kernel of a non-separable PPT state of rank four, acting on (C^ (D"^. 
The parameters p, q, r and s must necessarily be real. 

Proof. By dividing (1^ by 0i6 and (pS]) by 0Jg, which is possible according to 
016 7^ 0, we get 

(29) -^Vl + X^0'4 + -^-05 = = -|-V'l + :^V'4 + 

016 016 016 016 016 016 

Since {-01,-04, -05} is a linearly independent triple, the above equality implies Ai/0i6 = 
5i/0i6, A4/016 = ^4/016 and A5/016 = ■?5/0i6- In a similar way, from ([24]) and ([27]) 
we can get A2/026 = 6/026: -^4 / 026 = ?4 /026 and pA5 / 026 = P*C5/026> whereas 
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dSgi and ([Ml) give us A2/036 = ^2/(t>h, M/ft'ie = ^4/036 and q\z/(t)ze. = 9*6/036- 
From the equalities involving A4 and 6, we get 



(30) 



<^16 '^'26 </'36 

Together with X5l4>iQ — £,b/4>i%, the above equations give us A5/026 = ^5l4>2Q and 
A5/036 = ^/^f'Se- But 



(31) 



A5 ^5 . pAs P , ^ 

026 026 026 026 



In a similar way, from X^/cjizQ — S,zI4>^q and q\s,/4>m ~ 9*^5 / 036 '^an get q = 

q*. a 



3.5. An equivalence between generalized and orthonormal Unextendible 
Product Bases. In the following, we discuss item (4) of the list given in Section 
Let us start with a set of five vectors in (C^, 



(32) 



[ 



yil 912 (Pl3 914 <P15 
hi 022 023 024 025 
631 032 033 034 035 



and assume that any three of them are linearly independent, as in Proposition 13 . 1 3l 
For the moment, we do not require the vectors in p2p to be equal to 0i, . . . , 05 in 
(HU, but our ultimate goal is to apply the results we are going to obtain to ([T^. 
PSL (3, (D) transformations of the above set correspond to the multiplication of the 
matrix in (|32|) from the left by an element of SL (3, (D) and to the multiplication 
of the columns of ([32]) by arbitrary non-zero scalar factors. It is clear that we can 
PSL (3, (D) transform (|32|) to the following form. 



(33) 



^'l3 
^'23 
*33 



624 
6k 



ii5 
^'25 
^35 



By another PSL (3, (D) transformation, we get 



(34) 



1 
1 




1 



6'l3 

623 
633 



i'l4 



^25 
^'35 



^13 



^33 



f'24 
/'^,'4 



'^15 
6i'5 
6^'5 



We should remark that the matrix we multiply with from the left is well-defined, 
since 033 ^ according to the assumption about linear independence of triples. 
Let us transform once again, in the following way, 



(35) 



1 

1 
1 



Pl3 



^33 



ii'4 
i^'4 



■F13 



^33 



This is again possible because 035 ^ according to our assumptions. 



In a similar way as before, we see that 
fourth column by 1 



^24 7^ and 035 ^0. If wc multiply the 



24 and the fifth by 1/035, the above transforms to 



(36) 



1 a; ?; 
1 1 z 
< u 1 
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where we introduced the notation x := 



^13; 



24 1 



'34 /'/'24: 



'35. It is pretty straightforward to see that ah the coefficients x,y,z,t,u 



have to be different from zero according to the independent triples assumption 
Now, introduce the fohowing invariants 



(37) 



(38) 



Sl 



S2 




The numbers s\, S2 are indeed invariant. They do not change under the family of 
PSL (3, (D) transformations we were using in the consecutive steps (|32|) - (|36|) . Thus 
we can substitute 



(39) 



yii 

^•21 
^•31 



0\2 
^22 
^32 



yi3 

*23 
*33 



»14 
^•24 
^•34 



»15 
^25 
^35 



1 X 
1 
t 



in the above formulas for s\ and s^- In this way, we can easily calculate the values 
of the invariants, 



(40) 



Sl 



and 



S2 



UX 



Now, impose the conditions si > and > 0. From the ffist one, we clearly get 
u — —rz*, where r is a positive real number. Thus, we have the vectors 



(41) 



1 X 
1 
t 



y 
1 

-rz* 



Next, let us multiply from the left by a diagonal matrix diag ^1, 

as well as multiply the second column by 1 / \/\7r' , the fourth by 1 / 

fifth by ^/VP, where r' := rz* j z and stands for the square root of C G (D with 
the argument in [0, tt). Under such PSL (3, (D) transformation (|4ip changes to 




and the 



(42) 



1 x' 
1 
i' 



y 
1 



z 

1 



where z' is real and positive, and all the other parameters x' ,y' ,t' are non-zero. 
Moreover, the conditon S2 > transforms to 



(43) 



S2 



t'y' 



ux 



t'y' 

^ >0, 

x' 



simply by formula PO]) and the invariance of S2- The last equivalence holds by 
strict positivity of z' . In our next step, we we are going to multiply (j42p from the 
left by a diagonal matrix diag (Ci, C2, C3)i with CiiC2,C3 ^ and C1C2C3 — 1; and 
also multiply the consecutive columns, beginning with the first, by 1/Ci, I/C2, C4: 
(5 and Ce, where C4C5C6 7^ 0. Our aim is to choose the numbers Ci, . . . , Ce in such a 
way that ((42|l transforms to a set of vectors with orthogonality relations given by 
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a pentagon graph (that is, any two consecutive ones are orthogonal, and these are 
the only orthogonality relations). We would like to have 



(44) 



1 a b 
10 1a 
6 -a 1 



where a — z' and 6 is a positive real number in place of (|42|) . Let us write the 
numbers Q as rj-e*"^, where rj is a positive real number and aj E R. In order to 
obtain pi)) with a and b real and positive, certain phase matching conditions have 
to be fulfilled. Let us consider them first. If ay', at', ax' are such that y' = ry'e*"»' , 
t' — rt'C^"*' and x' — rx'C^"^' with Tx' , rt' and rx' real and positive, complex phases 
match correctly if and only if the following set of equations hold 



(45) 








a2 ^ 


- as 


= 


mod 27r 


(46) 








as ^ 


- ae 


= 


mod 27r 


(47) 








a2 ^ 


- ae 


= 


mod 27r 


(48) 








as ^ 


- "5 


= 


mod 27r 


(49) 






^5 + 


ai + 


a,y 


= 


mod 27r 


(50) 






"3 + 


a^ -f 


at' 


= 


mod 27r 


(51) 






Q!4 + 


«! + 




= 


mod 27r 


The 


requirement that C1C2C3 


— 1 adds 


a condition ai + 


equations (|45|l- 


gSl). 


However, 


a substitution of the form 


{cti 


,"2, as) ~^ 


(ai 


+ /3,a2 + 


5, as 


+/5) 


("4, 


"5,06) -> 



a2 + ^3 



mod 27r to 



(a4 - /3, as - /3, ae - /3) 

with an appropriately chosen /? can always bring ai + 0:2 + a^ to zero and it has no 
effect on ([151) - ([5T|) . Therefore, as long as existence of solutions is in question, we 
may neglect the additional condition. It is easy to see that the relations (|45)) - (l48l) 
are fulfilled if and only if q;2 = as = —as — — ag = —a mod 27r for some a G R. 
Thus the set of equations (j45 |) -([5T |l are reduced to 



(52) 



1 1 
-1 1 
1 1 



a 




-ay 


ai 




-at 






-ax 



mod 27r 



Interestingly, the 3x3 matrix in equation (|52|) has rank 2. A solution (a,ai,Q!4) 
exists if and only if 



(53) 



at' 



OLx 



mod 27r 



But this is exactly the positivity condition (j43p for the invariant si. Thus, if S2 > 
in addition to s\ > 0, we can cancel the complex phases, as in (j44|) . The only 
remaining thing to do is to match the modules, which gives us the following set of 
equations. 



(54) 



7'2?'5 



1, rare = 1, r2r6 = r^r^, r^riTx' = a, r^riry' = rar^rt' 



There is also an equation rir2r3 = 1, following from the requirement that C1C2C3 = 
1. As we see, there are five equations in ()54|) . and the variables ri, . . . ,re are six 
in number. Therefore, one can expect a solution to exist. It can easily be checked 
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that the following, with r e R, is a one-parameter family of solutions, 



(55) ri ^ . r, r2 = r, n = r, r4 = . , rs = -, = -. 

rx'ry' V r-^'Ty/ r r r 



By choosing r = 1/ ^arf /r^'ryi we can satisfy the additional condition rir2r3 = 1. 
Thus we have proved that the positivity of the invariants si, S2 guarantees that 
the family of five vectors ([5^ can be PSL (3, (D) transformed into ((44|) without 
permuting them. Obviously, a converse statement is also true, since the values of 
Si and S2 calculated from (|44|) are and b^/a^, respectively. In this way we arrive 
at the following 

Proposition 3.15. A set of five vectors {cxi}^^^ C (D^ with the property that 
any triple of them is linearly independent, can be PSL (3, (D) transformed, without 
permuting them, to the form (I44p with a and b real and positive, if and only if the 
invariants si and S2, defined in p7l) . are positive. 



Let us note that any set of five vectors {vi, . . . ,v^} C <D^ with orthogonality 
relations mods) — can be transformed by PSL (3, (D) transformations 

to the form (|33]). For sure they can be transformed to 



(56) [ Vi V2 W3 V4 V5 ] = 



1 a; y* 
1 1 a; 
y -X* 1 



with X and y complex. But since si = \x\^ > and S2 — \y/x\^ > in the above 
case, the argument following equation (j42|) tells us that a PSL (3, (D) transformation 
brings (1551) to the form (jH]). As a consequence, Proposition l3. 151 is a necessary and 
sufficient criterion for a set of five vectors 4>i,...,<j)z to be SL (3, (D) equivalent, 
without permuting them, to a set of vectors di , . . . , with orthogonality relations 

{Vi, Vi mod 5+l) = 0. 

From [2] we know that orthogonal UPBs in the 3x3 case always have five 
elements, and they are, up to permutations, precisely the sets of product vectors 
{vi (g) WiYi^^ ■^ith orthogonality relations {vi, Vi mod s+i) = and (^Wj,W(jj^i) mod 5+1) 
0. Consider the question whether an arbitrary set of five vectors {4>i ® V'i}i=i C 
C^ig) with linearly independent triples can be brought by PSL (3, C) x PSL (3, C) 
transformations to such {vi ® Wi\\^^, without permuting the vectors. In other 
words, what are the necessary and sufficient conditions for (pi ® ^i's to be con- 
vertible into Vi <Si Wi's with the orthogonality conditions given above. By using 
Proposition l3.151 we can already deal with the question about (j)i^s being convertible 
into Vi's. Namely, an PSL (3, (D) transformation on the first subsystem can bring 
the vectors {(t>i}i^i, without permuting them, to {fjjj^i with {vi,Vi mod 5-1-1) = 
if and only if the corresponding values of the invariants si and S2 are positive. 
We are only missing a similar criterion for ipi^s and w^'s. However, it is not dif- 
ficult to check that a permutation "^^^J 3 g 2 4^ brings any {wi}i=i 

with {wj,W(j + i) mod 5+1) = to = {^'^wlLl '^i*^^ (W;, mod 5+l) = 0' 

Therefore, it is sufficient to calculate the invariants ((37|) and ((38)) corresponding to 
the permuted vectors tpi ipaii) ^J^d check their positivity in order to tell whether 
the vectors ipi are convertible into some {wi}^^^ with the desired orthogonality 
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relations. Following the definitions ([ST]) and ([38]), let us introduce 
(57) S3 = 



and 

(58) S4 



1 ^/l 


V3 


oil 
W2 


1 ^/i 


05 V'4 














tp4 


• Ipl 

1 V'a(l) 


05 ^02 
V'<.(2) 


0^7(4) 


■ 1 V'a(l) 




V'<t(5) 








1 ^'-(l) 


V'a(2) 


V'o-(5) 


• 1 ^-(1) 




'0<t(4) 


1 V'l 


V'3 


■05 


■ V'3 


02 tp4 














V'2 


• V'3 
1 V'c.(l) 


05 "04 
V'.(2) 


0<t(3) 


'0<7(2) 


■0a(4) 


0,7(5) 








1 ^-(1) 


V'a(2) 


0<t(4) 


1pai2) 


V'a(3) 


0a(5) 



in accordance with [TT]. From the discussion above it follows that arbitrary five 
vectors tpi, . . . jijj^ in (D^ can be transformed, without permuting them, to {wi}^^^ 
with orthogonality relations {wj, W(^j_^i'f j^o^ 5+1) = if and only if the above invari- 
ants S3 and S4 are positive. Together with the previously obtained convertibility 
result between 0i , . . . , 05 and , . . . , U5 , the last result gives us the following 

Proposition 3.16. A set of product vectors {4>i <E) C C"^ ® (D'^ can be 

PSL(3,(D) X PSL (3, (D) transformed to an orthogonal UPB {vi®Wi}\^-^ with or- 
thogonality relations {vi, Vi mod 5+1) — and {wj, w^j^x) mod 5+1) = 0, without per- 
muting the 4>i (El ipi 's, if and only if the invariants si, S2, S3 and S4, defined in ([37]), 
(|38)) . (|57)) and (|58l) . are positive. 



Proof. Most of the proof has already been included above. Let {vi ® Wi}^^^ de- 
note an orthogonal UPB with the orthogonality relations {vi.Vi mod 5+1) = and 
(wj^w^j^i) mod 5+1) = for all i, j e {1, 2, 3, 4, 5}. The possibility to convert 



(59) 



■01 V'2 -03 04 V'5 



Vi V2 V3 V4 V5 



Wi W2 W3 Wi W5 



by PSL (3, (D) x PSL (3, (D) transformations, or by local equivalence in our usual 
terms, is the same as the possibility to separately convert {<t>i}^^i into {vi}^^i 
and {ipjj^^i into {wj}^^^ by some PSL (3, (D) transformations. However, we know 
that the first conversion is possible if and only if si and S2 are positive, while the 
second needs positivity of S3 and S4. Altogether, positivity of all the invariants Si, 
i = 1,2, 3, 4 is a necessary and sufficient criterion for the transformation ([59]) to be 
possible. □ 

In the context of product vectors in the kernel of a PPT state, as well as elements 
of an orthogonal UPB, permutations are obviously possible. Therefore we would 
like to have a version of Proposition 13 . 16l with no restriction on the ordering of the 
vectors {4>i (g) '0i}i'=i- 

Proposition 3.17. A set of product vectors {(f>i (E) C (D"^ ® (D'^ can be 

PSL (3, d) X PSL (3, d) transformed to an orthogonal UPB, if and only if for some 
permutation k the invariants si, S2, S3 and S4, calculated with the permuted vectors 
(t>K(i) and substituted for and ipi, respectively, are all positive. 
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Proof. Immediate given the fact [2] that an orthogonal UPB in a 3 x 3 system 
can always be brought by a permutation to a {vi (gi Wi}i=i with the orthogonality 
relations as in Proposition l3.16l □ 

Let us also note that, in accordance with [l^, not every single permutation of 
the five product vectors needs to be considered if we want to check whether they 
can be transformed into an orthogonal UPB or not. 

Remark 3.18. Only 12 permutations, given in Appendix \Al have to be checked in 
Proposition \3. 17\ in order to obtain a decisive answer. 

Proof. An explanation is included in lli and |13j . but we repeat it quickly here 
for completness. Let us denote by S5 the symmetric group of {1,2,. ..,5}. The 
permutations given in Appendix |X] are representatives of equivalence classes in 
55 of the regular pentagon subgroup G, generated by the cycle (12 34 5) and the 

inversion ^ ^ ^ ^ ^ 1 ^ ' '^^^ regular pentagon symmetry subgroup has the 

expected property that it does not change signs of si, S2, S3 and S4, just as it 
does not change orthogonality relations between the vectors {w;}^^]^ and {wj}^^^^. 
Therefore, we may divide S5 by G when we check positivity of the invariants in 
Proposition 13. 171 The number of invariance classes is 12 because #6*5 = 5! = 120 
and #G = 10. □ 

3.6. Determination of a PPT state by product vectors in its kerneL In 

the last part of the proof, we recall a number of surprising facts that were earlier 
reported in |131 Section 5] without a complete explanation. Here we fill in that 
little gap, and we collect a sufficient amount of knowledge to quickly obtain our 
main result. 

Note that, given a set of product vectors in Kerp, the conditions in Lemma 13.21 
are a set of linear equations for p. An idea, earlier presented in |13) . is to try to 
solve these equations assuming a specific form of the product vectors, namely (1141) . 
Let us repeat formula ([T4| here for the convenience of the reader. 

"10011 
1 1 p 
1 1 g 



(60) 



-01 1p2 i>4, '05 



10 11 
1 1 r 
11s 



We actually know from Proposition [STS] that there always exists a local SL (3, (D) ® 
SL (3, (D) transformation A (g) B that brings five vectors in the kernel of a non- 
separable PPT state of rank 4, possibly multiplied by some scalar factors, into the 
form (|60p with all triples linearly independent. Moreover, Proposition 13. 14l tells us 
that the parameters p, q, r and s are necessarily real numbers. By solving the linear 
conditions on a PPT state following from Lemma [3.21 with (pi ^ ipi, i = 1, 2, . . . , 5 
as in ((60)) substituted for <j) (g) ip, we will actually be solving a set of contraints 
on {A^^ (g) B^^)* p {A^^ ® ^^^)- However, according to the discussion in Section 
13.11 such local transformations are irrelevant to all the questions considered in 
this paper. Therefore we may simply assume that a PPT state p in question has 
the product vectors (|60|) in its kernel and check the consequences. As previously 
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reported by the authors of [13], the conditions {(pi ®^j,p {(pk ^i^i)) — Ofor i,j,k e 

) — imply the 



{1,2,3} together with p{(j) 
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04) = 
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with ai and bj real for aU i 
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, 6} and such that 




(62) ai+5i + 


62 = 


= 0, 


61 + 


0.2 


+ 63 
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03 + 64 
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(63) 63 + bi + 




= 0, 


65 + 


05 


+ 66 
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62 + be 
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(64) 










+ 62 


= 







Derivation of (|FT|) and the equations (|62I) - ([M)) is left as a simple exercise for the 
reader. It may be useful to consult Section 5.4 of ,13, in order to solve it. 

We still have not used the condition /9(05 <S) '4'5) = Oj which gives us additional 
six linear equations on ai, . . . , ag and 61, . . . , 65. 

(65) - r (61 + 62) + gr62 + sbi = 0, rbi - s {bi + 63) + psb^ = 

(66) —p (64 + 65) + 965 + psbi = 0, ^64 + 563 — ps (63 + 64) = 

(67) pb5-q{b5 + be) + qrbe^0, qbe + rb2 - qr {b2 + be) = 

Under the assumption of {(pi (g) V'i}^=i of the form (|60p being a gUPB, there exists, 
up to scaling by arbitrary real factors, exactly one solution to the equations (162^ - 
([67|). We know from Proposition l3.13l that the assumption is true for vectors (f>i^ipi 
in the kernel of a non-separable rank 4 PPT state in 3 x 3 systems. It is most 
important for us that there exist, up to scaling by arbitrary positive factors, exactly 
two solutions 



(68) ± 





qr-~s 
r(q-l) 
1 






r — s 
r{l-q) 







1 

s(l-p) 




r — s 
a{p-l) 











(r~s)(ps-q) 





p(l-s) 
r — s 
q-p 









s(p-l) 

r — s 
p(l-s) 



(p-s)(r-s) 
p{p-l)s{s-l) 











r — s 
q-p 





(gr~p)(r-a) 
q{l-q)(r-l) 



9(9-1) 














q{l-q)r{r-l) ^ 





r{l-q) 








q(q-i) 
(9-'-)('--s) 



The above matrix is well-defined since all the numbers p, q, r, s, p — 1, g — 1, r — 1, 
s — 1, p — q and r — s are nonzero as a consequence of all triples of vectors in (j60p 
being linearly independent. 

Note that, for both choices of sign, ([55]) is a symmetric matrix. Moreover, it 
is symmetric with respect to partial transpose. Therefore p is PPT iff it is positive 
definite. A necessary condition for (|68p to be positive definite is that all the nonzero 
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elements on its diagonal, as well as all nontrivial 2x2 minors of the form 
are positive. Altogether, we have six nonzero elements on the diagonal 



Pii 



Pij 
Pn 



(69) ± 



qr 



ps (r — s) {ps — q) {p — s) (r — s) 



^(9-1)' s{p-l)'p{p- 



q) {s-iy p{p-l)s{s-\y 

[qr — p){r — s) {r — q) {r 



q{p - q) [r - ly q{q - l)r [r - I) 



and six nontrivial minors 

(r — s) [qr ~ ps) 



(70) 



(9 ^ *) ^ s) '') ~ ^) 



r{p-\)s{q-\)' q{q-l)r{r -\y p{p 



{q - s)(r- s) 



(r — s) {qr — ps) 



1). (.-!)' 
(P -r){r- 



p{p-l){p-q)s{s-l)' p{p-q)q{r-l){s-l)' (p - q) q {q - 1) r {r - I) j 

The ± sign in (|69p corresponds to the choice we make in (|68p . We see that all the 
expressions in ([70| and (|69l) are quotients and products of the following nineteen 
numbers 

(71) p, q, r, s, p - 1, q - 1, r - 1, s - 1, p - q, r - s, 

(72) p — r, q ~ s, p — s, r — q, ps ~ q, qr ~ p, r — ps, qr — s, qr — ps 

Concerning the list ([7T|) . we already know that all its elements have to be nonzero. 
This follows from the condition of {0i ® '</'i}i=i being a gUPB. It turns out that the 
same holds for the elements of (|72p. The number qr ~ps must be nonzero, because 
otherwise the vector 

(73) (i}b®4>5- qr(l)4 ® -04 + g (s - r) 03 (g) ^3 - J' (p - 9) 02 '02 

would be of the form 0i (8) 0' + 0' <8) "01 , thus contradicting Lemma and Corollary 
13.121 In a similar way, one can show that p — r ^ and q — s 7^ 0. Let us now 
assume that ps — q = 0. In such case, we have the following submatrix in 



{r-s){ps-q} 
p(p-q){s-l) 







p(l-s) 
(p-s)(r-s) 
p(p-l)s(s-l) 



('^) =^ r-s (p-s)(r-s} -± 

p(l-s) p(p~l)s{s-l) _ 

In order for ([74|) to be positive definite for some choice of the sign ±, we need to 
have r — s = 0, which we know is impossible. Thus we have proved that ps — q ^ Q 
for p positive definite. Finally, the fact that qr ~ p, r — ps and qr — s must also 
be nonvanishing for p positive definite follows by a suitable modification of the 
above argument. Different submatrices need to be chosen, but otherwise the proof 
is identical. 

Our task in the following will be to relate positivity of all the numbers in ((69|) and 
([70)1 to the fact that all the invariants Si, . . . , S4, given in Section [XH are positive, 
possibly after we suitably permute the vectors 0^0 0^. Note that we already know 
that only the 12 permutations listed in Appendix [X] need to be considered. An 
explanation is included in the proof related to Remark 13. 181 Not to much surprise, 
the formulas for the invariants si, . . . , S4 for permuted vectors of the form (|60|) are 
always expressed as products and quotients including only the numbers listed in 
([TT]) . Explicit formulas can be found in Table [TJ To explain the notation we used 
in the table, it is sufficient to say, for example, that by using from Appendix lAl 
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0"! 
CT3 
0-5 

a^ 

0-11 



_P n — 1 ^-'i r 

1^ Q — p 1 — S 1 

g ' p ' s ' r~ 1 

\_ p—q 1 — r 1 

p ' g ' r ' s— 1 

p— 9 1 r s—r 

q ' q—1 ' s ' r— 1 

q—1 p ]_ 1 — s 

Q ' *?^P ' ^ ' r — 1 
q— p 1 _ s r-s 



0-2 
0-4 

(76 
(78 

0'12 



-£ n — 1 s 



s- 1,— 



1, 



p 

q 1— p r 

p—q ' Q—1 ' s—r ' 

q P-1 1 

1 — Q ' Q— P ' S— 1 ' 

P 1-g S 

Q— p ' p—1 ^ r — s ^ 



Table 1. Formulas for the invariants si, . . . , S4, calculated for vec- 
tors of the form ()60p permuted by the 12 incquivalent permutations 
of Appendix El 



to permute the product vectors (|60p . we obtain si = — p, S2 = (1 ^ -^s = r — 1 
and S4 = r/(s — r) as the expressions for the invariants. 

It turns out that the values of si,...,S4 corresponding to one of the permu- 
tations (7; have to be all positive to assure that p, given in ()68p . is a positive 
matrix for some choice of the sign ±. Our computer-aided proof of this fact 
consisted in simply checking all admissible sign choices for the numbers listed 
in ((7T|) and (|72p . We already know that neither of those numbers can be zero, 
and thus it seems that we have 2^^ cases to check. However, some further con- 
straints apply, which reduce this number considerably. First of all, the requirement 
that ±(p — s) (r — s)/{p {p — 1) s {s — 1)) of the list ([69]) and a very similar ele- 
ment (p — r) (r — s)/(p (p — 1) s (s — 1)) of (|70p have the same sign implies that 
p — r — ± (p — s), with the ± sign depending on the choice we made in (j68p . Along 
the same lines, by comparing the last element of (|69l) with the second element of 
(|70p . one can prove that r — q = =p(g — s). More importantly, the signs of the 
numbers listed in ((7T|) and (|72p are not all independent. Various relations have to 
hold between them. For example, p — 1 > clearly implies p > 0, and we cannot 
have a plus sign for p—1 and a minus sign for p. More sophisticated relations like 

(75) (gr - s > A p - 1 > A s < 0) => grr - ps > 

have to hold as well. We provide a more or less exhaustive list, consisting of 
76 elements, in Appendix [C] While some further relations could still possibly 
exist, the use of those listed in the appendix allowed us to confirm the necessity 
result mentioned above. When all the constraints are imposed, a comparably small 
number of 761 or 352 out of the 2^^ sign choices remain possible when or "— " 
is fixed in respectively. It then turns out that, by choosing an admissible sign 
configuration, all the numbers in the lists (|69p and ([70| can be made positive only 
if one of the quadruples listed in Table [T] consists solely of positive numbers. This is 
in full agreement with, and provides a rigorous, although not very insightful proof 
of the results reported in Section 5 of [T^. Actually, it turns out that there are 
precisely 12 admissible sign configurations that correspond to a positive p for some 
choice of the sign ± in ([55)1 and each of the quadruples in Table[T]is positive precisely 
for one of them. A complete list of the selected sign choices and the corresponding 
permutations is given in Appendix |B] Interestingly, 10 of them correspond to 
choosing the plus sign in ()68p . while only 2 to the minus sign. This is rather an 
uneven partitioning of the total of 12 configurations, which is somewhat puzzling. 
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To summarize, the computer-aided proof we carried out allows us to state the 
following. 

Proposition 3.19. A necessary and sufficient criterion for a generalized Unex- 
tendible Product Basis {(pi ^ ® belong to the kernel of a rank 4 

PPT state p is that there exists a permutation of the vectors (pi ® "04 th<^t H yields 
all the values of the invariants si, S2, S3 and S4, defined as in equations (|37p . (j38p . 
()57p and ()58p . positive. When checking positivity of Si, it is possible to consider 
only the 12 permutations, listed in Appendix [341 and the corresponding expressions 
for the invariants, given in Table [TJ 

Proof. First of all, let us note that a separable state p cannot have a gUPB in 
its kernel, since it must have a product state in its range. Thus in the following 
we may always assume that p is entangled. Let us prove sufficiency first. If the 
invariants are positive for the permuted vectors ^■^V'i '■= <P(y(i)®'4'a(i) : we know from 
Proposition l3. 171 that there exists a SL (3, (D)(8)SL (3, (D) transformation A^B such 
that the transformed vectors {A ^ B) (E) ip'^ = {A (E) B) <pcr(i) <S) 4'a(i) are elements 
of an orthogonal UPB {vi <Xi With no loss of generality, we may assume 

that the vectors Vi ® Wi are normalized to unity. In such case the projection 

5 

(76) p' := l-^lv, (^m,) (w,(8)w,| 

i=l 

has all the vectors Vi (E) Wi in its kernel and it is a PPT entangled state [J. The 
locally transformed PPT state p = {A (E) B) p' {A(E) B) has all the vectors (pi E) ipi 
in its kernel. 

In order to prove necessity, note that from the discussion above we know that pos- 
itivity of si, . . . , S4, possibly afer a permutation, is a necessary condition for a PPT 
entangled state p' with vectors (piE)tpi in its kernel to exist, provided that the vectors 
are as in equation ([60]). But any gUPB {(j)i (g) ■(/'i}i=i can be brought to the form ([60]) 
by a local transform, say C E> D. If we assume that a PPT state p has {pi (E 
in its kernel, then the locally transformed p" :— (C^^ E) -D^^)* p (C~^ (8) Z?"^) has 
(C E) D) (pi E) "ipi in its kernel. But (C E) D) (pi E) 4'i are of the form ((60| . From 
the above discussion, p" is PPT if and only if the invariants si, . . . , 54 are positive, 
possibly after we permute the vectors (C (E D) (j)i(Etpi. But C ®D does not change 
the value of the invariants, and thus <pi E) ipi, permuted in the same way as the 
(C ® D) (pi E) Ipi, must also have all of them positive. □ 

Let us also state the following result, which should be expected from the discus- 
sion above. 

Proposition 3.20. Let {(p^ (8) i/'ilLi C be a gUPB that yields, after a 

suitable permutation of the product vectors, positive values of all the invariants 
Si, . . . , S4. The PPT state p with {(pi E) *^ ^ts kernel is uniquely determined, 

up to scaling by a constant positive factor. 

Proof. We already know that the assertion of the proposition holds for gUPBs of the 
form We also know that any gUPB {(pi ® V'ili^i can be locally transformed 

so that it looks like in (pO]) . Let us denote the transformation which does it by 
C E) D. There cannot exist two PPT states pi and p2 with {(pi E) V'ili^i their 
kernels, because in such case the PPT states (C^i E) D-^)* pi {C-^ E) D-^) and 
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{C~^ (g) D-^y P2 {C-^ ® D-^) would both have the same gUPB of the form §0^ 
in their kernel, which we know is not possible. □ 

3.7. The main result. Using the knowledge from the previous sections, we can 
now easily prove our main result. 

Theorem 3.21. Positive-partial-transpose states of rank 4 m 3 x 3 systems are 
either separable or they are of the form 

(77) p= {AiS)B)* ^1 ~ ^ l'^' ® '^■') ® ^) 

with A, B E SL (3, (D) and {vi (E) Wi}^^^ an orthonormal Unextendible Product Ba- 
sis. In the latter case, they are entangled, and extreme in the set of PPT states. 
The rank of the partial transpose of the state is 4 in case of nonseparable states. 

Proof. In case of separable states, there is nothing to prove. Let p be a non- 
separable PPT state of rank 4 in a 3 x 3 system. We know from Proposition 
13.131 that there is a generalized UPB, say {4>i (S) ■'/'i}i=i: the kernel of p. From 
Proposition 13 . 1 9l we know that the corresponding values of the invariants si, . . . , S4 
must be all positive after we suitably permute the vectors Next, Proposition 

I3.16l tells us that there exists a SL (3, (D)(8iSL (3, (D) transformation A®B that brings 
{tpi (S) to an orthogonal UPB {vi (E) Wi}^^^. With no loss of generality, we 

may assume that the vectors Vi^Wi are normalized. From Proposition l3.20l we know 
that there exists, up to scaling, exactly one PPT state which has {vi ® ifili^i in 
its kernel. It must be 1 — Y^i=i ® '^i) (^j ®Wi\. The state given by the formula 
([77)) clearly is PPT, and it has all the vectors 4>i (8) ipi in its kernel. By using 
Proposition 13 . 201 again, we see that it must be equal to the p we started with. The 
fact that the rank of the partial transpose is 4 for non-separable states, is simply 
the assertion of Proposition [3791 □ 

In this way, we have obtained a full characterization of bound entangled states 
of minimal rank. Let us also mention a special property they have, which can be 
loosely described as saying that it is not enough for an entanglement witness to be 
indecomposable in order to detect them. 

Remark 3.22. According to Lemma 3], all PPT states of rank 4 m 3 x 3 
systems can be written as a sum of four projections onto vectors of Schmidt rank 
2. By Theorem \3.21[ or Provosition \3.9\. their partial transposes are also of rank 4 
and thus can be decomposed in an analogous way. Using the notation of \M3I, we 
can write that all such PPT states are elements of the cone S2.2- The dual cone 
822= T^2,2 consists of Jamiolkowski-Choi transforms of convex sums of 2-positive 
and 2-co-positive maps. Consequently, any entanglement witness that detects a 
PPT state of rank 4 m a 3 x 3 system is atomic f j7| /. This applies in particular 
to the witness discussed in Example 1 of U8\j . 

3.8. Another use of product vectors. From the long argument above we have 
certainly learned that product vectors in the kernel play an important role for PPT 
states. It seems desirable to give an additional example of how they can be used 
to explain properties of PPT states reported in ^4| . We have the following 

Proposition 3.23. A PPT state p acting on (D^ (E) C" must have a product vector 
in its range. 
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Proof. Assume that there are no product vectors in R{p). According to Lemma 
I2.13[ which originates from [16; Lemmas 1 & 2], the kernel of p, which is R{p)'^, 
contains a set of product states {(l>i (g) V'i}i=i such that their partial conjugates 
</>* (gi ipi span (D^ igi CD". But this is impossible in our case, because according to 
Lemma [3.11 all the vectors (j>* ^ ipi must be in Kerp-^i. If they spanned (D^ ® CD", 
we would necessarily have p^^ = as a consequence of positivity of the partial 
transpose. Thus we see that R (p) must contain a product vector, or in other 
words, it is not a CES. □ 

One can check in Table I of [14^ that all PPT states numerically studied by 
the authors in the 2x4 and 2x5 cases actually have a product vector in their 
range. The above proposition shows this is a general phenomenon. On the contrary, 
Tables II - IV of [14 reveal PPT states, supportes on 3 x n and 4x4 subspaces, 
that do not admit a product vector in their range. 

4. Conclusion 

It seems profitable to sometimes drop the orthogonality condition for Unex- 
tendible Product Bases, just as it is useful to extend the group of allowed trans- 
formations for PPT states from local special unitary to local special linear trans- 
formation jllll25) . By using this approach, with help of Bezout's Theorem from 
algebraic geometry, we have obtained a characterization of bound entangled PPT 
states of minimal rank, discovered earlier in numerical searches by the authors 
of [11]. The states are all SL (3, (D) ® SL (3, (D) transforms of projections onto or- 
thogonal complements of orthogonal Unextendible Product Bases. In this sense, 
they are a "deformed" version of the bound entangled states introduced in [T] . Cer- 
tain methods we used apply to more general cases, and we hope they can help to 
explain even more peculiarities of the tables in |14) . We already did explain a few 
of them, including the absence of non-separable PPT states p of rank 4 with the 
rank of p^^ different from 4 and the fact that 2 x n PPT states must always have 
a product vector in their range. The use of product vectors in the kernel and the 
range of a PPT state is a common trait to most of the work reported above. Note 
that the paper [30j by other authors discusses extreme PPT states in the 2x4 case, 
which can be regarded as another step towards an improved understanding of the 
results of [Tl] . 

One of our main aims was to show that algebraic geometry can be useful in 
solving very concrete problems in quantum information theory. We could probably 
say that we went somewhat further than the papers |20 |[22l[23] . but in a similar 
spirit. Apart from solving the main example in Section |3l we showed, for instance, 
that some of the findings of [T6j can be found in a disguised form in a textbook 
like [15]. We also used the opportunity to collect basic results about general Un- 
extendible Product Bases, which seem to be missing from present literature. The 
gUPB framework arises in a natural way when local SL transformations are used 
instead of local unitary ones, and we could see this in Section |3| 
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Appendix A. Inequivalent permutations 



This appendix gives a list of representatives of the 12 equivalence classes of 
the symmertic group S5 under left multiplication by the regular pentagram group, 



generated by the cycle (12 3 4 5) and the axial symmetry 

refer to them in Remark 13 . 1 81 and in the discussion in Section 
compliance with |13) . but we provide it here for completness 



0-1 



(T4 



^2 : 



CT5 : 
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(Til : 
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erg 



1712 



. We 



3 4 5 
3 2 1^ 
The list is in full 



Appendix B. Sign choices that yield a positive p 

In Table [21 we provide a complete list of twelve admissible sign choices for the 
numbers listed in equations (I7T|) and (17^ that yield a positive matrix (|68p . It turns 
out that for each sign configuration in Table [21 there is exactly one permutation ai 
among those listed in Appendix [X] that yields a positive value of all the invariants 
Si, . . . , S4 defined by formulas (|37l) . (|38| . (|57l) and (|58| (for more details, cf. Section 
I3.6|) . It should be mentioned that the rows in Table [2l corresponding to cg and cio 
need a minus sign in ()68|) . while all the other ones correspond to choosing a plus 
sign. To explain the notation we used in the header of Table [H it is probably 
sufficient to say that the symbols pp, qq, rr and ss denote 1, <?— 1, r — 1 and 
s — 1, respectively, while pq, rs, pr, rq, qs, qrp, qrs, psq, rps and qrps stand for 
p — q, r — s, p — r, r — q, qr — p, qr — s, ps — q, r — ps and qr — ps, respectively. 



Appendix C. Sign constraints 

In this appendix, we give a list of constraints on the numbers listed in equations 
(|71]) and (|72|) . They are contained in Tables [3] and [H Some of the constraints are 
simple relations like l>0^p>0, but there are also more sophisticated ones. 
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Table 2. Sign choices that yield a positive p and obey ah the 
constraints of Appendix [C] 



For example, 

(78) (r<OAg-l<OAr-ps>0)^(/r-ps>0. 
Alternatively, the above formula can be written as 

(79) ^ {r < A q - 1 < Ar - ps > A qr - ps < 0) . 



This is also the convention we adopt in the following tables. For example, formula 
((79| corresponds to the following row in Table |3l 
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An explanation of the notation we used in the headers of the tables can be found 
in Appendix [Bl 
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Table 3. Non-admissible 



sign choices. Part I. 
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Table 4. Non-admissible sign choices. Part II. 



